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Abstract 

Using the fact that the energy eigenstates of the equilateral triangle infinite well (or billiard) are 
available in closed form, we examine the connections between the energy eigenvalue spectrum and 
the classical closed paths in this geometry, using both periodic orbit theory and the short-term 
semi-classical behavior of wave packets. We also discuss wave packet revivals and show that there 
are exact revivals, for all wave packets, at times given by T^ev = 9fia'^/AhTr where a and fi are 
the length of one side and the mass of the point particle respectively. We find additional cases 
of exact revivals with shorter revival times for zero-momentum wave packets initially located at 
special symmetry points inside the billiard. Finally, we discuss simple variations on the equilateral 
(60° — 60° —60°) triangle, such as the half equilateral (30° — 60° — 90°) triangle and other 'foldings', 
which have related energy spectra and revival structures. 
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I. INTRODUCTION 



The study of the connections between the quantized energy eigenvalues of a bound state 
and the classical motions of the corresponding classical point particle has undergone some- 
thing of a renaissance as the ability to experimentally probe the quantum-classical interface 
has dramatically improved. Theoretical methods such as periodic orbit theory 3|; 0, 
for example, provide very direct connections between the energy spectrum and the closed 
trajectories of the classical system. The time- dependence of fully quantum mechanical wave 
packet solutions of the Schrodinger equation also depends critically on the energy spectrum, 
both for the short-term semi-classical propagation as well as for longer-term, purely quantum 

mechanical effects such as wave packet revivals. Such revival phenomena have been observed 

n 

in a wide variety of physical systems, especially in Rydberg atoms [^], and calculations exist 
for many other systems j^. 

Two-dimensional billiard systems have provided easily visualizible examples relevant for 
)oth types of analyses. For example, the periodic orbit theory analysis of square/rectangular 
^ and spherical/circular geometries were among the first performed, while on the ex- 
perimental side, the energy level structure and statistics of microwave cavities ^ (modeling 
2D billiards of arbitrary shape) have been probed to test both periodic orbit theory and to 
obtain statistical evidence of chaotic behavior. Measurements of conductance fluctuations in 
ballistic microstructures have been tentatively used to identify frequency features in the 
power spectrum with specific closed orbits in a circular (and stadium) billiard. More recently, 
the realization of atom-optics billiards with ultra-cold atoms in arbitrary shaped 2D 
boundaries confined by optical dipole potentials, has allowed the study of various chaotic and 
integrable shapes such as the stadium, ellipse and circle, again for short-term, semi-classical 
propagation. 

The one- ^| and two-dimensional square wells, with their integral energy 

spectra, have provided simple examples of exact wave packet revivals, wherein initially 
localized states which have a short-term, quasi-classical time evolution, can then spread 
significantly over several orbits, only to reform later in the form of a quantum revival in which 
the spreading reverses itself, the wave packet relocalizes, and the semi-classical periodicity is 
once again evident. Circular billiards have also been recently studied ^| where approximate 
revivals have been found to be present under some circumstances. 
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One billiard system which has not been analyzed in as much detail is the case of the 
equilateral (or 60° — 60° — 60°) triangle quantum infinite well or billiard. If one thinks of 
the square billiard as an = 4 regular polygon, with the circle being the N ^ oo limit, 
the = 3 case of an equilateral triangular billiard is an obvious extension. Since the exact 
bound state wavefunctions for the corresponding infinite well problem have been available 
for some time ^|, ^3]; our purpose here will be to apply these known results and 
analyze the quantum mechanics of this system in the context of periodic orbit theory as well 
as wave packet revivals. Given how few quantum systems, especially billiard geometries, 
have closed-form solutions, it is very useful to have another such example analyzed in as 
much detail as possible. In addition, the energy eigenvalues for this case, described by 
simple quadratic powers of quantum numbers, just as for the 2D square billiard, will be seen 
to yield exact quantum revivals, providing another interesting 'testbed' for possible future 
experimental tests. 

In Section II, we review the solutions of the equilateral triangular well, wavefunctions and 
energy eigenvalues, and make some comments about simple variations resulting from various 
'foldings' of the standard equilateral triangle billiard. We then perform a periodic orbit 
theory analysis of this system in Sec. Ill and discuss the exact wave packet revivals present 
in this system in Sec. IV, as well as visualizing the short-term semiclassical propagation 
wave packets and their relationship to classical closed orbits. 

II. ENERGY EIGENFUNCTIONS AND EIGENVALUES 
A. Results for the equilateral triangle billiard 

The energy eigenvalues and position-space wave functions for a particle of mass fi in an 
equilateral triangular infinite well (or billiard) of side a have been derived, in a variety of 
different contexts, by at least three different groups 



3|i -^o^ definiteness, we will 
assume a triangular billiard with vertices located at (0, 0), (a/2, y/3a/2), and (—a/2, v^a/2), 
as shown in Fig. 1. The resulting energy spectrum is given by 

£'(m, n) = [m'^ + n^-mn) (1) 

for integral values of m, n, with the restriction that m > 2n. (In what follows, we will use 
the notations of Ref. [l6[ for the energies and wavefunctions: we use n for the particle mass 



to avoid confusion with various quantum numbers.) For the case of m > 2n, there are two 
degenerate states with different symmetry properties ||l6| which can be written in the forms 

/ 2Txny\ 



and 



^(Lln)(^'l/) 
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Extending earher resuhs, we have here included the correct normahzations, since we will 
eventually expand Gaussian wave packets in such eigenstates. 

For the special case of m = 2n there is a single non-degenerate state for each n, given by 



^(2i,n)(^'l/) 
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Clearly these states satisfy 
',{±) 



^(Z',n)i-^^y) = ±^(t]n)i^^y) 



(o) 



and the ip^^^2n n) y) states also satisfy 



{m=2n,n) ^ 

',{±) 



0. 



The authors of Ref. [ig] note that the V^l^ {x, y) states also satisfy such relations as 



(±) 



and 



(4) 



(5) 



(6) 



(7) 



and several other similar ones. The extremely symmetric character of the (2n, n) states can 
also be seen by the fact that it can be rewritten in the form 



i'^$n,n)i^^y) 



8v/2 f2'Kny\ ( -nniy — \/2>x)\ ( 7m(y + \/3x) 



(8) 



These wavefunctions satisfy the Schrodinger equation with the energy eigenvalues in 
Eqn. as well as vanishing on the boundaries defined by the lines y = ±^/3x and 



to their x —x parity properties. The authors of both Ref. |15| and [16| argue that the 
solutions above also form a complete set of states by explicitly calculating the Weyl area 
rule for the energy level density 



and showing that this relation is saturated by the energy spectrum in Eqn. (0). While we 
will focus on the case of the equilateral triangle billiard, it is useful to note that several 
simple variations on this system (see below) can also be derived from these results. 

We note that these solutions, while derived in the context of the explicit triangular 
area shown in Fig. 1 (shown in bold), are also solutions, when extended by reflections, 
in the other triangular regions (flipped in sign) and hence over the complete hexagonal 
area including reflections in the horizontal axis. (This connection, in fact, was used in the 

n 

explicit construction of the solutions in one approach [1^.) Thus, these solutions (with 
their integrally quantized energies, will also form a subset of the solutions for the hexagonal 
(A^ = 6 regular polygonal) billiard. 

B. Comparison to the square and circular billiards 

The degeneracy pattern of energy eigenvalues for the equilateral triangle billiard is very 
similar to that for the square (A^ = 4 polygon) and circular (A^ = oo polygon) systems. For 
the square billiard (of side a), the energy eigenvalues are given by 



y = \^a/2. The eigenfunctions corresponding to different eigenvalues can be shown to 
be orthogonal, while the two with the same eigenvalues (^[^^ Jx, y)) are obviously so due 





(9) 




(10) 



with eigenfunctions given by 




(11) 



with 
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and any integral m,n > 1 allowed. For m = n, there is a single state, while for m ^ n there 
is a two-fold degeneracy. The analogs of the degenerate (±) states of the equilateral triangle 
billiard for the square infinite well are given by 

^(m!n)(^'2/) = ^ [Un{x)Um{y) + (m ^ u) (13) 

^(nln)(^'2/) = Un{x)Un{y). 

This form is useful since it allows one to discuss the energy eigenvalues and eigenfunctions of 
the 45° — 45° — 90° triangle [loj billiard formed by 'folding' the square along a diagonal , 
since the '4^\jnn)i^->y) satisfy the appropriate boundary condition along the new hypotenuse. 
Additional foldings along diagonals are also possible [3| and the energy spectrum can be 
used to analyze these cases as well. 

A similar pattern of degeneracies, also useful for analyzing 'half-well' problems, is found 
in the circular well [2^. In this case, the wavefunctions in polar coordinates can be 
written in the form 

V;(r, 9) = [NrJim\{kr)] [Nge'""^] where m = 0, ±1, ±2, ... (14) 

where the radial wavefunctions are normalized via 

/ [NrJim\{kr)Y rdr = l. (15) 
Jo 

The states corresponding to non-negative values of +m and — m are degenerate. A standard 
set of angular wavefunctions can be written as 

\Nee'"''^f de = 1 or A^^ = (16) 







but for the purposes of visualizing the solutions of bound state problems, real wavefunctions 
are also possible, given by 

cos{m9) / ^/ti m > 1 
ei^j(^) = <J sin(m0)/v/? m > 1 . (17) 
l/V^Ti m = 

Thus, there is a doubly degenerate set of states labeled by (m > 1, n,.) for m > 1, and a single 
set of states (all of the s-state wavefunctions) with (m = 0,77.^). One of the {m > l,nj.) 



sets can then be used to describe the energies of the 'half-circular' billiard, as either the 
cos{m6) or sin(m6') solutions (or some linear combination thereof) can be used to satisfy the 
boundary condition along an infinite wall on any diameter, giving a 'folded half- well'. 

C. Variants on the equilateral triangle 

Using the explicit wavefunctions for equilateral triangle, it is easy to construct solutions 
of the 'half-triangle' (namely the 30° — 60° — 90° triangle) obtained by folding the equilateral 
triangle along one bisector: the odd-parity solutions in Eqn. for example, explicitly 
satisfy the new boundary condition along the bisector given by x = 0. The resulting energies 
are still given by Eqn. (^, but now with only a single copy of the m > 2n states allowed, 
which agrees with earlier analyses The wavefunctions in the new 'half- well' are 

identical to the V'(„^„(a;, y) except with an additional normalization factor of y/2 to account 
for the geometrical 'footprint' of the smaller billiard. 

The three lowest-lying such states, corresponding to (m, n) quantum numbers and en- 
ergy values of (3,1), (4,1), and (5,2) and 7Eq,13Eq, and 19Eq respectively (where Eq = 
(/i^/2/ia^)(47r/3)^) are shown on the top line of Fig. 2. 

It is easy to shown that if one lets (m',n') = {2m, 2n) then all of the energy eigenvalues 
are multiplied by the same factor of 4 and that there are nodal lines which split up each 
'half-triangle' into four equal parts corresponding to a simple folding. The corresponding 
states for this subset of quantum numbers are also shown in Fig. 2 (in the bottom row). 
This possibility is part of a general pattern where the mapping (m', n') = {fm, fn) yields 

'copies' of the (m, n) wavefunctions (for both the even and odd states) contained inside 
the original triangular billiard; one can see this more explicitly by noting that 

^(/mjn) i^^ = ^(S„) (2^' 2/; «//) • (18) 

Slightly less trivially, the quantum number mapping given by (m', n') = (2m — n,m — 2n) 
gives states which are a factor of 3 times in energy, with nodal lines along the perpendicular 
bisectors, as shown in (the middle row of) Fig. 2. A second application of this map gives 
{m" ,n") = (3m, 3n) and a resulting common increase by a factor of 9, as mentioned above. 
Both sets of results can be then be easily understood in terms of repeated foldings of the 
'half-triangle' onto itself, with increasingly small areas and larger energy eigenvalues scaling 



7 



in the appropriate manner, just as in Ref. for diagonal foldings of the square billiard. 
For this case, we might note that the odd wavefunctions explicitly satisfy the relation 

V^{2m-n,m-2n)(3^' «) = «/3) • (19) 

Other simple linear transformations of quantum numbers of the form (m', n') = {am + 
(3n, 7m + Sn) also return allowed quantum numbers in the triangular well, corresponding to 
higher energy states, but without the same obvious geometrical significance. A few examples 
of such mappings are discussed briefly in the Appendix. 



III. PERIODIC ORBIT THEORY ANALYSIS 



A. General background 

As mentioned above, the notion of density of energy levels has already been applied to the 
solutions of triangular billiard to the extent that the averaged energy level density, Po{E), in 
Eqn. ©, has been evaluated to confirm completeness, as in Refs. H and Q- It is perhaps 
then an obvious next step to expand on this type of approach by using a periodic orbit 
theory analysis 

Q. a. 0. 

In this context, the energy level density is naturally split into 
the smooth, slowly varying part {po{E)) and a second, oscillatory term which is dominated 
(in a saddle point approximation of the path integral) by classical orbits whose actions, 
S^{E), correspond to periodic orbits or closed paths. Specifically, using the quantized energy 
eigenvalues, labeled collectively as En, one writes 



00 00 r / s" (E) \ 

S(E - En) = p{E) = poiE) + Y,Y1 P^^P P ( - ) 

n = ^ •n= 1 'V I- ^ ' 



(20) 



n=l p=l 7 

where each periodic orbit is characterized by a label 7 and p = 1,2, ... runs over all pos- 
sible repetitions of such trajectories. For billiard systems, this expression can be Fourier 
transformed and simplified in terms of 

»+oo 00 



^ f + OO ^ 

p{L) = J2 '^") ^'^^ dk = J2' 

n=l -f-"^ n=l 



to give 



p{L)=Y.Y.P^^vKL-pL,). (22) 

p=i 7 
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Thus, an evaluation of p{L), using the quantized fc„ values given by the En = h^kl/2p, will 
exhibit a series of 5-function like 'spikes' at multiples of the lengths of the 'primitive' closed 
paths, i.e., at L = pL^. For numerical studies, an evaluation of 

N 

p^{L) = Y,e''-' (23) 

n=l 

for large, but finite will give a series of sharp peaks corresponding to the closed orbits, 
with the smooth po{E) piece corresponding to a large, but uninteresting feature at L = 0. 



B. Periodic orbit theory analysis of the equilateral triangle billiard 

For a numerical evaluation of pn{L), we require the quantized wavenumbers obtained 
from the energy spectrum in Eqn. (^, namely 

k{m,n)a = m? + 11? - mn (24) 

so that 

{k{m,n)] m/2 > n > 1} = {7.255, 11.082, 11.082, 14.510, 15.102, 15.102, ...} , (25) 

where including the correct degeneracy, as we will see, is important. An analysis of pn{L), 
evaluated using Eqn. (|^. then requires knowledge of the possible closed paths or periodic 
orbits in the triangular billiard. 

A particularly easy way to visualize the possible classes of closed orbits in this geometry 
is by making use of the fact that the 2D plane can be tiled with a triangular lattice, as 
illustrated in Fig. 3. Consider a specific equilateral triangle billiard, in this case the one in 
the lower left hand corner. By repeated foldings, the entire 2D plane can be covered, and 
'equivalent' triangles are shown with appropriate points identified. Any point in the original 
triangle can then be connected to the corresponding point in another 'identified' partner (as 
illustrated by the arrows) and the resulting path inside a specific triangular billiard can be 
obtained (visualized) by repeatedly folding the two triangles until t hey overlap. (We note 
that this simple construction has been used for the square billiard [2l| to obtain the path 
lengths for the closed orbits.) In Fig. 4 we illustrate several cases corresponding to different 
starting points within the triangle, yielding identical overall path lengths, so long as the 
initial angles are the same. We note that, just as for the square and circular billiards, there 



are an infinite number of possible closed orbits of the same path length, here corresponding 
to different initial locations of the orbit. For example, the trajectories shown in Figs. 4(a) 
and (b) both correspond to path lengths of L = 3a (and multiples thereof), while those in 
Figs. 4(d), (e), and (f) all yield path lengths of ^/Sa. Fig. 4(c) shows a unique (isolated) 
orbit which, because it closes on itself in a very symmetric way, corresponds to a single 
possible closed orbit with path lengths given by multiples of 3a/2. 

The general expression for the allowed closed orbit lengths can then be obtained using 
this construction. For example, without loss of generality, we can find the path lengths 
between a single corner in the initial triangle in Fig. 3 (say the corner at (0, 0), labeled by 
the circle) and any other corresponding corner. For those points corresponding to triangles 
which have been folded into 'even' rows, we can, for example, use the points 

Xij^ = ii{3a) and Uji = ji{y^ci) (26) 

which gives 

d{ii,ji) = ^J{x,,-0y + {yj,-0y 

= a^9i? + 3jH^V9(2ii)2 + 3(2i0^. (27) 

This corresponds to initial angles (with respect to the horizontal) of 

tan(^) = !^ = (28) 

Xii llVo 

For points in triangles which have been folded into 'odd' rows, we similarly have 

= (2i2 + l)y and y^.^ = (2j2 + l)^ (29) 
giving closed path distances 



2 

and initial angles 



d{i2,j2) = V(^*2-0)2 + (%,-0)2 

" V9(2i2 + 1)2 + 3(2j2 + 1)2 (30) 



tan(^) = ^ = J^^i±l) . (31) 
Xi, (2^2 + l)\/3 



We can combine these two cases into a single expression, namely 



d{ij)^^\/9(ir + m' (32) 
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where are either both even or both odd. The corresponding angles are given by 

tan(0) = ^ . (33) 

The path length features for these classes of closed orbits (and their recurrences) which 
have L = d{i, j) < 20a, along with the corresponding [i, j) labels and angles, are shown in 
Table I. This construction finds all of the allowed closed orbits, except for the special case 
described by Fig. 4(c). 

We can rewrite Eqn. (jH^ in a more symmetric form by defining 

{p + q) = i, ip-q)=J or p= , q = . (34) 

Since i,j are either both even or both odd, the p,q defined in this way are integer valued 
and we find that the path lengths can be written as 



d{p,q) = aV3\/p'^ + pq + q"^ (35) 

and this form will be useful when we discuss how the closed or periodic orbits arise in the 
context of wave packet propagation in the next Section. 

We then evaluate Pn{L) (using the = 1000 lowest lying eigenvalues) and plot |pAr(L)p 
versus < L/a < 20 in Fig. 5. The path length features collected in Table I are shown 
as vertical dotted lines and the vertical scale has been chosen to emphasize these features. 
We note also that there are two much smaller features, indicated by arrows, corresponding 
to the isolated, single closed orbit trajectories at the first two odd multiples of L/a = 3/2, 
namely L/a = 1.5 (denoted by (i, j) = (2, 0)') and L/a = 4.5 (for 3(2, 0)'). Even recurrences 
are, of course, hidden under the standard = (2,0)) L/a = 3.0 orbits, while higher 

odd recurrences are close to other features corresponding to different orbits: for example the 
5(2,0)' feature at L/a = 7.5 is near a more prominent (5, 1) feature at L/a = 7.55, while 
the L/a = 10.5, 7(2, 0)' feature is shadowed by the L/a = 10.53 (7, 1) orbit. 



C. Periodic orbit theory analysis of the 'half-triangle' 

We can easily extend the analysis of the equilateral triangle billiard to the 'half-triangle' or 
30°— 60° — 90° triangle by using the reduced set of (m > 2n) energy eigenvalues corresponding 
to the 'ip\J„\{x,y) states which vanish along the y-axis as discussed in Section II C and the 

I I ft/ ^ IL J ^ ' 
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evaluation of Pn{L) proceeds just as above. The only substantiative changes in the allowed 
classical closed orbits can be visualized as has been done for various foldings of the square and 
circular wells in Ref. Q|. For all allowed orbits (save one special one), the introduction of 
an infinite wall along a bisector simply serves to reflect all trajectories back onto themselves 
inside the half-triangle, giving identical overall path lengths for closed orbits, as exemplified 
in Figs. 6(a) and (b) for two of the (i, j) = (1, 1) orbits. In this case only, there is a second 
isolated path (shown in Fig. 6(c)) with half the expected path length, namely L/a = \/?>a/2 
(and multiples thereof) which we will call the (1, 1)' orbit. 

We thus expect to see in the \pn{L)\'^ spectrum the same path length features as for 
the full-triangle, but with much smaller new features at odd multiples oi L/a = 
or L/a ^ 0.866,2.598,4.330,6.062,7.794,9.526,... and so forth. We plot, in Fig. 7, pn{L) 
for both the full- well (dashed) and half- well (solid) cases, adjusting the overall scales to 
emphasize the similar shape, and we note the agreement at all of the standard, non-isolated 
path length features in Table I, as well as for the isolated L/a = 3/2 (2,0)' orbit. For the 
half-well data, however, we also see clear evidence for all of the new isolated (1, 1)' orbit 
features (shown by arrows.) 



IV. TIME-DEPENDENCE OF GAUSSIAN WAVE PACKETS 
A. General time-dependence of wave packet solutions 

The connections between the energy eigenvalue spectrum of a quantum system and the 
classical dynamics of particles can also be explored through the time-dependence of wave 
packet solutions of the Schrodinger equation. For a localized wave packet in a system 
characterized by a single quantum number, n, an arbitrary wave packet solution can be 
written in the form 

oo 

V^(x,t) = ^a„V^„(x)e-^^(")*/^ (36) 

n=l 

One then typically expands the energy eigenvalues (assuming integral values) about the 
central value, hq, used in the construction of the wave packet to write 

E{n) ^ E{no) + E'{n^){n - no) + \E"{no){n - n^f + ^E"'{no){n - n^f + ■■■ (37) 
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in terms of which the classical period, revival, and superrevival times 
tively by 

27rh „ 2TTh „ 2'Kh 



ll\ are given respec- 



\E'ino)\ 



T 

J. ri 



\E"(m)\/2 



super 



For systems with two quantum numbers 



3, 



\E"'{no)\/Q 



(38) 



and energies labeled by E{ni, 71,2), the 



patterns of both classical periodicities and revival times can be much richer as we will see 
below. We first examine how the closed orbits for the equilateral triangle in Eqn. ()32|) arise 
in this formalism. 



B. Short-term, semi-classical periodicity in the equilateral triangle billiard 



For a system, such as the 2D billiard, for which the quantized energies depend on two 



quantum numbers, (m, n), there are two corresponding classical periods 
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T, 



(m) 



27rh 



T. 



in) 



27ih 



I given by 

dE{rn,n)/dm ' dE{m,n)/dn' ^^^^ 

Classical periodic orbits will arise when the ratio of these two times is a rational number, 
or equivalently when 



{po) 



d -H^ci (40) 
which defines the classical period, t'^"\ for such closed orbits. For the equilateral triangle 
billiard, with energies given by 



E[m, n) 



2/ia2 



47r 

T 



(m^ + — mn) 



we find that 



pim) 
-cl 



9fia 



21 



{2m — n) 



pin) 
-cl 



{2n 



m] 



The condition leading to closed orbits in Eqn. PH]) can then be written as 



(2m — n) 



in) 



m] 



(m) 



P 



or 



n 



m 



2p + q 
2q + p 



(41) 



(42) 



(43) 



If we substitute this condition into the energy spectrum in Eqn. (j4ip . as well as equating 
the quantum energies with the classical kinetic energy, ^Vq/2, in the billiard (where vo is 
the classical speed) we are led to the association 



E{m, n) 



"~9~ 



2/ia^ 



2 , 2 / 2p + g 

m + m 



2q + p 



m 



2p + q 
2q + p 



(44) 
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or 



This implies that 
'2/it>oa 



m 



(2g + p) 



m 



3(p^ +pq + q^) 
(2g + p)2 

/ 2fiVQa 



n 



{2p + g) 



The period for classical, closed/periodic orbits is then given by Eqn. (jl0|l as 



(45) 



(46) 



(po) 



(m) _ aySyp^ +pg + _ L{p, q) 



d 



where 



(47) 



(48) 



L{p,q) = d{p,q) = a\/'i\/p^ +pq + q^ 

and we have used the alternative form of the path lengths for periodic orbits in Eqn. ()35j] . 
(A similar analysis for the circular billiard can be found in Ref . [l^ , making use of the WKB 
approximation for the non-integral energies found there, while the corresponding derivation 
for the square billiard is entirely straightforward.) 

To visualize this result, we make use of a specific form for Gaussian wave packets, namely 



where 



i){x,y;t = 0) = ipo{x;xo,Pox,b)ipo{y;yo,Poy,b) 
'ilJo{x;xo,pox,b) 



1 



^ipOx(x-xo)/h^-(x-xo)'^/2b^ 



with a similar expression for ipodi; yo,Poy, b). 

The initial expectation values for the x variables are given by 

6' . b 



(49) 
(50) 



and 



{Px)o= POx 



4 + 



{pDo = pIx + 



2 

262 



V2 



n 



V2b 



with similar results for y. The expectation value of total energy is 



-{Pl+Pl) 



{POxf + {POyf + 



62 



(51) 



(52) 



(53) 



2m ^ 2m 

Just as with the allowed wave functions, wave packets in the billiard must satisfy the 
appropriate boundary conditions at the edges of the well. However, the simple Gaussian 
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form above can be used to describe an initial state so long as the initial location, {xo,yo), 
is well away from the edges of the potential well, so that a negligibly (exponentially) small 
error is made. In that case, such a Gaussian form can be easily and reproducibly expanded 
in terms of eigenstates, namely 

oo oo 

ij{x, y-t = o) = («(2n)^{i„)(^' y) + «{i)^(™In)(^' y)) 

m>2n n=l 

oo 

+ E«(ln)Cr.)(^'^)- (54) 
n=l 

For a sharply peaked Gaussian, we can (very accurately) approximate the required expansion 
coefficients by extending the required integration of (x, y) over the triangular region to one 
covering all space, making use of the exact form of the wavef unctions. Given the simple 
results in Eqns. (0), (0) and (jH), we can perform the necessary Gaussian integrals in closed 
form using 



Cx\ , 6V2^ 
a 



j^^-iCxo/a^-b-^{-C/a+po/hf/2 



and 



g.po(x-xo)Ae-(-xo)V2ft^ sinf— 1 dx = ^ 

a J 2i 



^iCxo/a^-bHC/a+po/h)y2 ^gg^ 



_^-iCxo/a^-b^-C/a+po/h)y2 



We can then construct wave packets with arbitrary initial locations (xo,yo) inside the well 
(again, so long as they are well away from any edge, typically more than a few times 
Axo = Ayo = b/V2) with arbitrary initial momenta, {poxiPoy), which we parameterize via 
(po cos(6'),po sin(6')) to make connections with the specific angles characterizing closed orbits 
in Eqn. (jH^ and Table I. 

All quantum wave packets have a characteristic spreading time, so in order to exhibit 
behavior typical of a semi-classical regime, we want to have classical periodicities shorter 
than a typical time over which a Gaussian wave packet will spread significantly. For the 
free-particle form in Eqn. ()5()|1 . one can write 

Axt = Axo^/TTJt/W where to = ^ = ^^{Axo)\ (57) 
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The classical periods for closed orbits will be given by Eqn. ()47|) and the two shortest such 
periods will be 

'ijnin) = — = — (58) 

and we want to have Tj:[°\min) < to- For purposes of illustration then, we have taken the 
following numerical parameters: 

h,2fx,a = l , b = 1/1072 so that Axq = Ayo = 0.05 , and po = 1500. (59) 

With these values, to = (0.05)^ = 2.5 x 10-^ while T^[°\min) = (0.58,1.00) x lO^^ and 
one can see many classical periods for most of the lowest path length classical closed orbits 
before significant spreading occurs. 

In order to probe whether a wave packet has returned to something close to its initial 
state after some expected classical period, it is standard to examine the autocorrelation 
function i22l], defined, for example, in one-dimension as 

A{t)= I ^*{x;t)^{x;0) = y"\an\^e'^-'/^ (60) 



/ + 00 
n=l 



which also explicitly depends on the energy spectrum. Evaluating A{t) for our two- 
dimensional wave packets, we then plot, in Fig. 8, \A{t)\ versus t/r (where r = cl/vq) 
for Gaussian wave packets which are initially located at the geometric center of the triangu- 
lar well ((xo,?/o) = (0, v^fl/3)), and with the initial momentum in the 9 direction: we then 
show results for various 6 values in the range 0° < < 30°. (Because of the symmetry of 
the billiard, such plots will be symmetric about 6 = 30°.) 

We expect to see features in A{t) at values of 9 and classical periods L{p,q)/vo given 
in Table I as we 'scan' up in angle. The stars in Fig. 8 denote the location (in {6,Tj^°'') 
parameter space) where we would expect classical particles to return to their starting point 
in phase space and we do indeed see well-localized structures in A{t) near each expected 
location. (The dotted curves which connect the 'stars' representing the exact closed orbits are 
simply there to guide the eye.) We see this same pattern if we initially localize the Gaussian 
wave packet to places in the well other than at the geometrical center, corresponding to the 
infinite sets of non-isolated classical orbits which all have the same path length. 

One exception to this occurs for parameters near the initial location (xq, yo) = (0, ^/Sa/A) 
and with 6 = and the results for A{t) for such this case are shown (on the bottom line) 
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as the dotted peaks. They coincide with the standard case, but also exhibit structure at 
classical path lengths given by odd multiples of 3a/2 and are due to the special isolated 
= (2.0)' path shown in Fig. 4(c). 
Such wave packet solutions can also be constructed in the 'half-well' case, using the ap- 
propriately normalized versions of the V'(m^„)(3^) v)- As mentioned above, an initial Gaussian 
wave packet will be expandable (to exponential accuracy) in this way, provided it is contained 
within an appropriate 'fiducial' area of the triangle, well away from any of the boundaries. 
We note that because such a wave packet would also be expandable in the original 'full-well' 
(with coefficients given by Eqn. (jMl)), the expansion coefficients in the 'half-well', ^l^''^), 
will necessarily satisfy the relation 

aS") , = V2aS"^ , (61) 

(m,n) (m,n) V / 

Since the must saturate the complete probability for the 'half-well' version of the wave 

packet, we must then have 

oo oo oo oo 

m>2n n=l m>2n n=l 

and the odd expansion coefficients (for such wave packets which are well inside one half 
of the original 'full well') give exactly half of the total probability. For wave packets in 
the 'half- weir, the only new closed orbit path length feature, beyond those seen in Fig. 8, 
arises from wave packets which are constructed to mimic the special isolated (1, 1)' orbit in 
Fig. 6(c). 



C. Long-term revival behavior 

The long-term behavior of wave packets can include the possibility of revivals, wherein 
initially localized states, after undergoing a number of semi-classical periods eventually 
spread out during a so-called collapsed phase, only to reform later into something very 
much like the initial state. The time scales governing these revivals, in a two-quantum 



number system are given by 



0, Q 



{d^E/dm^)/2 ' {d^E/dn^)/2 ' (d^E/dmdn) ' ^ ^ 
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For the equilateral triangle under consideration here, we have the simple result 

nn{m) _ nn{n) _ nn{m,n) ^ nn _ 9/^^ (f>A\ 
■^rev ~ -^rev ~ rev — -^rev — y'J'^J 

and exact revivals (just as for the 2D square well) are present in this system, with a single 
revival time guaranteed for any and all possible initial wave packets. 

To visualize this phenomenon, we again make use of Gaussian wave packets, with the same 
parameters used above, but with vanishing momentum {pq = 0) where there would be no 
obvious classical periodicity. We locate several such wave packets at initial positions (keeping 
away from the edges) along the y-axis, as shown in Fig. 9, and plot the autocorrelation 
function over a time equal to one revival time, Trev In each case we do find an exact revival 
at the expected time, but we also notice special situations with exact revivals at shorter 
times, due to additional symmetries of the eigenfunctions, which result in the vanishing of 
certain expansion coefficients. For example, for (xo,?/o) = (0, v^a/3) (geometric center of 
the well), one finds revivals at multiples of Trev/^ since only expansion coefficients for which 
(m + ra) is a multiple of 3 are non-vanishing, due to the nodal pattern of the eigenstates 
superimposed on the circular symmetry of the {pox,Poy) = (0,0) Gaussian wavefunction; in 
a similar fashion, if {xo,yo) = (0,v^/4) (initial position half way down the bisector), one 
finds exact revivals at multiples of Trev/4:. Even at these special points, if the momentum is 
made non- vanishing, only the true revivals at Trev remain intact, since the purely radially 
symmetric nature of the wavefunctions is lost. (Similar cases such as these are also seen in 
the ID and 2D square well jl^ for zero-momentum wave packets when obvious symmetries 
dictate that certain subsets of expansion coefficients vanish.) 

Finally, since the 'half- well' version of the triangular well also has integrally related 
energies, exact revivals are also present in the 30° — 60° — 90° triangle, but without the 
special cases mentioned above present. 



V. CONCLUDING REMARKS 



We have used the triangular billiard, with its simple energy spectrum, as a model system 
to probe periodic orbit theory and wave packet revivals. Similarly to the cases of the 
square {N = 4 polygon) and circular (A^ — * oo polygon) billiards, we find that for the 
equilateral triangular (A^ = 3) infinite well the periodic orbit theory analysis agrees well 
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with expectations, including special isolated orbits for various 'folded' versions; for the case 
here, we find a special isolated orbit in the 'full-well' version as well. The short-term, 
semi-classical propagation of quantum wave packets shows evidence, in the autocorrelation 
function, of the same closed orbits, including the isolated ones. 

The equilateral triangle billiard also exhibits exact wave packet revivals (due to the 
integral nature of the energy eigenvalues), just like the 2D square billiard, as does the 'half- 
well' version, for all wave packets. 

The case considered here might also be used as a basis for a more systematic examination 
of the quantum energy spectrum of the = 6 polygon, namely the hexagonal billiard. The 
solutions here will form a subset of the allowed spectrum of the = 6 case and an exami- 
nation of both the periodic orbits and wave packet revivals in such a system will require the 
complete enumeration of the energy eigenvalues, either in closed form, or perhaps using var- 



ious approximation methods 



23|: having a subset of energy eigenvalues and wavefunctions 



for which the results are known exactly would, of course, provide a very useful 'benchmark' 
for any numerical calculation. That case would also be another interesting example of an/N 
billiard {2^. Such a result would be useful in probing the pattern of wave packet revivals as 
one goes from the = 3, 4 cases (where exact revivals are present) to the N 00 (circular) 
case, where only approximate revivals exist. 
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APPENDIX A 



In Sec. n(c) we exhibited several simple transformations of the (m, n) quantum num- 
bers of the equilateral triangular well which, when applied to allowed (m, n) pairs, returned 
quantum numbers and corresponding energy eigenvalues appropriate to higher energy states 
which, in turn, could be easily visualized in terms of simple 'foldings' of the half-well. For 
example, the maps (m', n') = [fm, fn) gave energies which are times larger, correspond- 
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ing to p smaller 'copies' of the initial {m,n) state inside the full triangular well. The 
less obvious mapping given by (m', n') — (2m — n,m — 2n) gave energies which were three 
times larger, corresponding to obvious 'foldings' in thirds, as evident from the resulting wave 
function plots, as in Fig. 2. 

If we define a dimensionless function, e(m,n), via 

E{m^ n) = Eoe{m, n) = Eo{m^ — mn + n^) (Al) 

where Eq = (^^/2/ia^)(47r/3)^, we can see that 

e{fm, fn) — pe{m, n) and e(2m — n,m — 2n) — 3e(m, n) . (A2) 

This second mapping is a special case of a general set of such mappings which return subsets 
of the allowed energy eigenvalue space. If we pick any pair of integers, (p, q) , of the same 
form (p > 2q) as the allowed quantum numbers themselves, and define the transformation 

(m', n') = 7(^,5) [m, n\ = {pm - qn, {p - q)m - pn) (A3) 

one can easily confirm that 

e(m', n') — [p^ — pq-\- q^] {vn? — mn + n^) — e{p, q)e{m, n) (A4) 

and we find all such energies in the spectrum of the equilateral triangle. Two consecutive 
applications of the same transformation then give 

(m", n") = T(j,^q) [T^p,q)[m, n]] = ((p^ - pq + q^)m, (p^ -pq + q^)n) (A5) 

so that 

e{m'\n") = [e(p, g)]^e(m, n) . (A6) 

One can similarly show that two successive transformations with different labels, in either 
order, give rise to equivalent sets of quantum numbers in the sense that the corresponding 
energies are identical, namely if 

(m", n") = T(p,g) [T^r,s) [m, n\\ or (m", n") = r(^,,) [T(p,g) [m, n\\ (A7) 

both give 

e(m", n") = [e(p, q)e{r, s)]e(m, n) (A8) 
20 



which also appear in the original energy spectrum. We note that in order to see that 
the (m", n") given in the two orderings are, in fact, identical, one sometimes must use 
the {m,n) {m,m — n) equivalence in Eqn. (|7j). These transformations do not seem to 
have the same obvious geometrical significance as the T(2,i)(m, n) three way 'folding' or the 
(m',n') = {fm,fn) 'copying'. 
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Tables 



(^J) 


e (deg) 


path lengths [Lja < 20) 


(^,J) 


d (deg) 


path lengths {Lja < 20) 


(2,0) 


0.0 


3.00, 6.00, 9.00, 12.00, 15.00, 18.00 


(11,5) 


14.7 


17.06 


(13,1) 


2.5 


19.52 


(4,2) 


16.1 


6.24,12.49,18.73 


(11,1) 


3.0 


16.52 


(9,5) 


17.8 


14.18 


(9,1) 


3.7 


13.53 


(5,3) 


19.1 


7.94,15.87 


(7,1) 


4.7 


10.53 


(11,7) 


20.2 


17.58 


(12,2) 


5.5 


18.08 


(6,4) 


21.1 


9.64,19.28 


(5,1) 


6.6 


7.55,15.10 


(7,5) 


22.4 


11.36 


(13,3) 


7.6 


19.67 


(8,6) 


23.4 


13.08 


(8,2) 


8.2 


12.12 


(9,7) 


24.2 


14.80 


(11,3) 


8.9 


16.70 


(10,8) 


24.8 


16.5 


(3,1) 


10.9 


4.58,9.16,13.75,18.33 


(11,9) 


25.3 


18.24 


(13,5) 


12.5 


19.97 


(12,10) 


25.7 


19.97 


(10,4) 


13.0 


15.39 


(1,1) 


30.0 


1.73, 3.46, 5.20, 6.93, 8.66, 10.40, 
12.12,13.86,15.59,17.32,19.05 


(7,3) 


14.0 


10.82 



Table I. Path lengths (from Eqn. (j32|) ) and initial angles (from Eqn. (j33|) ) for closed orbits 

(and recurrences) with L/a < 20. 
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Figure Captions 



Fig. 1. Geometry of the equilateral triangle billiard (shown in bold) with vertices at (0,0), 
(+a/2,y3a/2), and (-a/2, y3a/2). 

Fig. 2. Nodal patterns for the lowest-lying energy eigenstates for the 'half triangular' well 
(using the u) i^i Eqn. Q) on the top row. (The white (black) areas correspond 

to positive (negative) values of the wavefunction.) The bottom row corresponds to 
{m',n') = {2m, 2n) states which are a factor of 4 larger in energy, while the middle 
row shows those related via (m', n') = (2m — n,m — 2n) which are a factor of 3 higher 
in energy, illustrating various 'foldings' of the half-well. 

Fig. 3. Geometrical visualization leading to the construction of the path lengths corresponding 
to various closed orbits in the triangular billiard. 

Fig. 4. Closed orbits corresponding to several = (1,1) cases (bottom row) with path 

lengths given by multiples of -\/3a and for (i, j) = (2, 0) cases (top row) with path 
lengths given by multiples of 3a, except for (c) which is an isolated orbit (described 
as (2, 0)' in the text) with path length given by 3a/2 (and multiples thereof.) 

Fig. 5. Plot of |piv(-^)P versus L/a for the triangular billiard, evaluated using the 1000 lowest- 
lying energy eigenvalues. All of the expected path length features from Eqn. and 
Table I are indicated by vertical dotted lines. The small bumps indicated by arrows 
correspond to the isolated (2, 0)' orbits in Fig. 4(c) with path lengths given by multiples 
of 3a/2. 

Fig. 6. Closed orbits in the 'half-triangle' well for the = (1, 1) case. For (a) and (b), the 
path lengths are still given by L(l, 1) = \/3a (and multiples thereof), but the 'folding' 
induces a single isolated (1, 1)' orbit with L = \^a/2 as shown in (c). 

Fig. 7. Plot of \pn{L)\'^ versus L/a, for the triangular billiard (dashed curve) and the 'half- 
triangle' billiard (solid curve). Vertical dashed lines indicate the expected path length 
features in Table I, while the vertical arrows indicate the small new features (at L/a = 
V3/2 and multiples thereof) present in the 'half-triangle' due to the new, isolated 
closed orbit in Fig. 6(c). 
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Fig. 8. Autocorrelation function, versus t/r, (where r = a/vo) for Gaussian wave 

packets with the parameters discussed in Eqn. Wavepackets initially localized 

at (xo,?/o) corresponding to the geometric center of the triangle are given momenta 
(po cos(0),Po sin(0)) and \A{t)\ plots are shown for 9 values in the range (0°,30°). The 
'stars' indicate the location of classical closed orbits (corresponding to the values in 
Table I). The dotted \A{t)\ data for the 6' = 0° case corresponds to the initial wave 
packet localized at {xq.uq) = (0, v^a/4) which corresponds classically to the isolated 
(2,0)' orbit in Fig. 4(c). 

Fig. 9. Autocorrelation function, |A(t)P versus t, for zero-momentum Gaussian wave packets 
initialized at various points along x = in the triangular well. There are exact 
revivals for all cases shown (where |y4(Tret,)| = 1), but also special, shorter-period 
revivals at multiples of Trev/^ and Trev/4: for wave packets initially at positions of 
special symmetry in the well, namely at the geometric center (0, ^/3a/3) (square box) 
and at (0, V^a/A) (half way down the bisector.) 
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